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As is well known, the fluctuations from a stable stationary nonequi- 
librium state are described by a linearized nonhomogeneous Boltzmann- 
Langevin equation. The stationary state itself may be described by a non- 
"^ ' linear Boltzmann equation. The ways of its linearization sometimes seem 

^ I to be not unique. We argue that there is actually a unique way to obtain 

J^ ■ a linear equation for the fluctuations. In the present paper we treat as an 

O I example an analytical theory of nonequilibrium shot noise in a diffusive con- 

ductor under the space charge limited regime. Our approach is compared 
with that of Schomerus, Mishchenko and Beenakker [Phys. Rev. B 60, 
C^ . 5839 (1999)]. We find some difference between the present theory and the 

approach of their paper and discuss a possible origin of the difference. We 
^_^ . believe that it is related to the fundamentals of the theory of fluctuation 

^5 ' phenomena in a nonequilibrium electron gas. 

o 

> 

^ '. I. INTRODUCTION 

(N 

Q . The present paper is devoted to the theory of shot noise in space-charge limited dif- 

fusive conduction regime. The motivation of the present paper can be formulated as 
follows. It is well known that the fluctuations from a stable stationary nonequilibrium 
state are described by a linearized nonhomogeneous Boltzmann-Langevin equation (see, 

i-rt ■ for instance, [|1]-|31)- At the same time, the stationary state itself is described by a non- 

linear Boltzmann equation. There are instances where the ways of linearization of the 
nonlinear Boltzmann equation seem to be not unique. We believe, however, that in each 
such case there is a unique way to obtain the linearized Boltzmann equation for the fluc- 
tuations and we give general considerations to find out this way and indicate it for the 
particular case treated in the present paper. 

We work out a theory of nonequilibrium shot noise in a nondegenerate diffusive con- 
ductor under space charge limited regime. This regime is extensively discussed in the 
literature (see, for instance, Refs. [§3). The current noise under such a regime has been 
studied recently by Monte Carlo simulation by Gonzalez et al. |jlO|. Quite recently the 



noise under the same conditions has been studied analytically by Schomerus, Mishchenko 



and Beenakker |TT[]. Their general finding was that due to the Coulomb correlation be- 
tween the electrons the shot noise is reduced below the classical Poisson value. Both the 
authors of Ref. [|1^] and of Ref. [|ll| came to the conclusion that under certain conditions 
the suppression factor in the nondegenerate 3D case can be close to 1/3. 

Later on Nagaev [jl2[ has shown on a special example that unlike the 1/3 noise reduc- 
tion in degenerate systems, the noise suppression by the Coulomb interaction is nonuni- 
versal in nondegenerate systems. The noise suppression in such systems may depend on 
the details of electron scattering. 



We agree with the conclusion |lTO| , pJ| that the reduction of the shot noise power in 
nondegenerate diffusive conductors can be sometimes close to the value 1/3 predicted 
theoretically for a three-dimensional (3D) degenerate electron gas. As we have mentioned 
above, we also come to some conclusions that may prove important for the general theory 
of fluctuations in nonequilibrium systems. As is well known, the fluctuation phenomena 
in nonequilibrium stable systems are described by a linearized Boltzmann equation. We 
would like to use the example analyzed in detail in the present paper to show that one 
should be careful performing the linearization. In particular, there is a difference between 



the analytical procedures used in Ref. [jTT| and in the present paper for calculation of 
the shot noise power. We discuss the origin of such difference and its implications. As 
the point leading to the discrepancy is very subtle it demands a rather detailed analysis 
that we perform in the present paper partly repeating the calculations of Ref. |jll| with 
some modifications. Our starting point is the Boltzmann equation that is formulated for 
description of the stationary state and will be applied for analysis of the fluctuations. 



II. BOLTZMANN EQUATIONS 



We will consider the simplest model, exploited in Ref. W^, for the diffusion-controlled 
and space charge limited transport. As the starting point we use the Boltzmann equation 
in the presence of electric field 



Qi+^Pjfp 



0, 



(2.1) 



^p/p - (v|: + eE A + jp j /p. (2.2) 

Here we have introduced the collision integral Jp describing the electron scattering: 

Ipfp = E i^P'pfp - WpP'fp') . (2-3) 

p' 



(we deal with the nondegenerate statistics, so that /p ^ 1). 

Sphtting the distribution function into the even and odd parts with respect to p one 
gets 

Jp ~ 7^ l/p ^ J-p) ■ 

We assume that the colhsion operator acting on the even (odd) part of the distribution 
function gives an even (odd) function. This may be due either to the central symmetry 
of the crystal itself and the scatterers or to the possibility to use the Born approximation 
for calculation of the scattering probability. The first split equation is 

df- (9f+ (9f+ 

l + ^lF + '-^l^-Wp-. (2.4) 

Being interested in relatively small frequencies of fluctuations UTp ^ 1 where Tp is the 
characteristic value of I~^ we can neglect the time derivative and express /_ as 

/; - -^;' (vf + eEvg) . ,2,5) 

Inserting this expression into the second split equation for /+ ^ f{e,r,t) and averaging 
over the surface of constant energy in the quasimomentum space we arrive at 



f .df ( d 



dt \ dxr 



+ eE^^^j iy{e)D^,{e) (^ A + eE,^^ / = - E ^(^ " ^r^)!^^ (2-6) 



where the term on the right-hand side describes the inelastic collisions while the density 
of states z/(e) and the diffusion tensor Dapie) are defined as 

u{e)D^p{e) = Y.6{6-ep)vJ-'vp, u{e) = Y.6{e - Sp). (2.7) 

p p 

The electric field obeys the Poisson equation 

roo 

KV-E = 47re[n(r,t) -n"'], n(r,t) = / de iy{e)f{e,r,t), (2.8) 

Jo 

where k, is the dielectric susceptibility, n'^'^ is the equilibrium concentration (equal to the 
concentration of donors). In what follows we neglect it as compared with the nonequilib- 
rium concentration n. 

The part of the distribution function contributing to the current consists of two terms 
proportional to the spatial and energy derivatives of f{e,r,t) respectively 

Jc = eE v/p = -eu{e)D^f,{e) (^^ + eE^-^^ f. (2.9) 



Let us consider the case Dts ^ L^, where L is the sample length, r^ is the energy 
relaxation time (of the order of J^™*^^) ). Then we can omit the term on the right-hand 
side of Eq. ( |2.6|) which describes the energy relaxation. Under the same conditions we 
get the following Boltzmann equation for the fluctuations of the distribution function (we 
remind that we consider here low frequency fluctuations with uo <^ I-p k. l/r-p where t^ 
is the characteristic time of elastic collisions) 

+ eE^— 5jZ + e5EZ^j^ = ey^e, x) (2.10) 



\dxr, "del •'" "9£ 

where 



^J!^ = ej^^^a^/p = gZ- ep{e)D^p{e) 



d ^ d 

+ eEp- 



Sf^ + eSE^^A (2.11) 



dxf3 de 

and where the source of the current fluctuations g'^ is related to the Langevin forces y^ 
as 

9Z = eY.5{e~e^)vJ^X^ (2.12) 

p 

y^(£,x)=5^5(£-£p)2/^ = 0. (2.13) 

p 

The last equality is a consequence of elasticity of scattering that leads to the particle 
conservation within the surface of constant energy in the quasimomentum space. 
The correlation function of the Langevin forces is well known [^ 

< yp(r)yp.(r') >^= (Jp + Jp05rr'<^pp'/p. (2.14) 

Integrating Eq.( |2.1CI| ) over e we obtain the continuity equation 

A— / de5i^{e,x) = —6J^{x) = 0, (2.15) 

dx JO dx 

which implies that the low frequency current fluctuations are spatially homogeneous. 

III. DISTRIBUTION FUNCTION. 

Consider a semiconductor with a uniform cross section A connecting two identical 
metallic electrodes. The sample's length L is assumed to be much bigger than the elastic 
scattering length / and much smaller than the inelastic one. We use the ID versions of 
the Boltzmann equations describing evolution of the distribution function along the dc 
current direction. 



To obtain the stationary solution of Eq. (|2.6|) in the accepted approximation we rewrite 
it in the form 

(^ + "^^)^'^"'^)=^^^)^'^"^- ^^-^^ 

Here we assume that the current density at x = 0, j{£), does not vanish only for e > 0. 
In the absence of tunneling j{e) at the contact x = should have the property 

j{e) ^0 for T ^ 0, (3.2) 

T being the temperature. This condition should be valid, irrespective to whether one has 
a Schottky barrier or an Ohmic contact. The total current J given by Eq. ( |3.3| ) below 
should have of course the same property. 

The solution of Eq. ( |3.1| ) is a function of the total energy S 

S = e + U{x), 

where U{x) = eip{x) — eip{0). It can be found employing, for instance, the inverse differ- 
ential operator 

Ox Jo 

One has 

"'*'^'^'* " d +eE(x)d ^'^^^■^^^^ " exp[eip{x)ds]-^exp[-eip{x)ds]6{x)j{e) = j{S) 

and j{e,x) has at a given x nonzero values only if e > —U{x), {£ > 0). The total 
current through the sample is 

/■OO poo POO 

J = A dej{e,x) = A de j[e + U{x)] = A dSj{S). (3.3) 

Jo J-U{x) Jo 

Now, from Eq. (|2.9|) we get 

or 

/(e,.) = -,|fl/>j3^^4^^ + /[£] (3.5) 

where X{e) = i'{e)D{e). We have taken into account the boundary condition at the 
source. Equation (|3.5|) can be rewritten as 



f[£-U{x),x] 



f[s- u{L)] r d^/x[s - f/(0] + f{s) r d^/X[S - f/(0] 

Jo Jx 



\^/x[s-u{0] 



(3.6) 



where j{e) is expressed through the difference of the distribution function at x = and 
X = L 



An advantage of the form we have chosen for Eq. ( p.6| ) is its physical transparency. 
The first term on the right-hand side gives the contribution from the right boundary 
while the second one gives the contribution from the left boundary. The solution clearly 
demonstrates that the thermally exited carriers injected from the contact at x = L make 
negligible contribution to the distribution function f[S — U{x),x] since f{S) ^ f[S — 
U{L)] {£ > 0), for the parameter \U{L)\/k-QT is assumed to be large. Neglecting this 
term in our solution of Eq. ( |3.6| ) we come to the solution already obtained in |^ by 
assuming absorbing boundary conditions at the current drain. 

IV. FIELD DISTRIBUTION 

We use the Poisson equation to determine the selfconsistent electric field that can be 
expressed through the obtained distribution function. We consider such values of x that 
X > X-, where 

-U{xe) > ^ ~ A;bT 



^ei 



roo 

/ d£ v\E - U{x)\ 
Jo 



/•oo rL 

xf[£-U{x),x]= dSu[£-U{x)]j{£) 



d^ 



-urn 

Finally we get 

K I d^U J fL d^ 



' ' '''' 'Jx eX[£-U{0] 



47re2 iy[-U{x)] dx^ eA Jx X [-f/(0] 



(4.2) 



Let us check that for large x this equation is consistent with the requirement of a 
uniform total current. Assuming u(e) = uqe'^''^^^ and D[e) = Dqc'^^^ we integrate Eq. ( |2.9|) 
over the transverse coordinates and energy 



-e— / dei^{e)D{e)f{e,x) + ^^ ^-[-Uix)Y—E^{x) 4.3 

dX Jo lOTT dx 



We integrate by parts the second term and take into account that at x > x^ we can 
neglect S as compared to |f/(a;)| and use the Poisson equation ( p.Sl) . The first term in 
Eq.( [4.3|) can be simphfied in the same way [note that due to Eq. (|3.5| ) the distribution 
function f{e, x) has nonzero values only for e > —U{x) ] 



/OO /"OO 

de u{e)D{e)f{e, x) = d£ u[£ - U{x)]D[£ - U{x)]f{£ - U{x),x) 
-U{x) Jo 

= D\-U{x)] dSu\S-U{x)]f{S-U{x),x)=D\-U{x)] -E (4.4) 

JO 47re dx 



Here in the second equality we have taken into account that S ^ |f/(a;)|. Inserting 
Eq. ( [4.4|) into Eq. ( [4.3|) we get the following simplified equation 



2 



It can be used to verify the self-consistency of our approach. Indeed, multiplying Eq. (O) 



by f/*+'^/2 and taking the derivative we arrive at Eq. (|4.5| ) that has been obtained from 



the equation for the current. A dimensionless version of Eq. (4.5) is 



X' {^XY - XX'") = 1 (4.6) 



where the dimensionless potential x is related to ip by 



^-i^m^r^^^"^^- 



V. CURRENT AND FIELD FLUCTUATIONS 

In what follows we consider the particular cases s = 0, D{e) = Dqe; s = —1/2, 
D(e) = Dqe^^'^; and s = 1/2, D{e) = DqE^^'^. We begin with investigation of the energy- 
independent-scattering-time case s = 0. This case can be related to the scattering of 
electrons by the neutral impurities, such as hydrogen-like shallow donor and acceptor 
states. The scattering is analogous to the scattering of electron by a hydrogen atom [l^ 



(with the effective Bohr radius a^). The scattering cross-section turns out to be about 
27rh/{paB) times larger than the geometrical cross-section vra^ (that would result in an 
energy- independent scattering time). 

In the case of defects with deep energy levels we encounter a short-range scattering 
potential with the scattering length about atomic length. The scattering cross-section 



does not depend on the energy. As a result, the scattering rate is proportional to the 
electron density of states e^/^ and the diffusion coefficient v'^t is proportional to e:^^^, i.e. 
s = — f /2. (This is one of the main scattering mechanisms in metals since the scattering 
length there is determined by the radius of screening which is of the order of interatomic 
distance.) The case s = —1/2 (that, in particular, describes elastic scattering by acoustic 
phonons) and s = 1/2 will be discussed at the end of this section. 

A. Energy independent scattering time 



Integrating Eq. ( p.ll| ) over e, we get 



4 i^J^ - G.) = -e^ r deu{e)D{e)6Ue,x) + ^^-^E{x)6E^{x) (5.1) 

A ax J-u(x) oTT ax 

Note that the Fourier transform of the current fluctuations 5J^ due to Eq.( |2.1CI| ) is spatially 
homogeneous. Here G^ is the current fluctuations source integrated over the energy and 
transverse coordinates 

G^{x)= dedr±g^{£,r) (5.2) 

Jo 

< G{x)G{x') >^= eM de de' ^ 5{e - ep)5{e' - ^pO^^^x' 

■^0 JQ pp, 

11/- 

X — -— / dr^drf± < y^y^, >^ (5.3) 



'P -'P 



The odd (with respect to p ^ — p) part of the distribution function vanishes after one 
inserts it into the correlation function ( 2.14| ) of the Langevin forces and subsequent inte- 



gration over p and p'. As a result, we are left with the integral of the even function 

< G{x)G{x') >^= 4,. < G\x) >,, (5.4) 

{G\x))^ = 2e^A / defie, x) J2 <^(^ " ep)v,—v, = 2e^A / de v{e)D{e)f{e, x). (5.5) 
Jo ^ Ip Jo 

The second term on the right hand side of Eq. ( |5.1| ) can be simplifled in the same way as 
Eq. O 

/ deu{e)D{e)5Ue,x) = D{-Uix))-——SE^ (5.6) 

JO 47re dx 

and flnally we get the equation for 6E^ 



ax \ ax I 2 ax 



An 



ADqk 



[oJuj — Guj). 



(5.7) 



In order to justify the simplification in Eq. ( |5.6| ) we will show that 5fc^{e,x) is also a 
function acquiring nonzero values only at e > —U{x). Indeed, from Eq. ( |2.10| ) and 
Eq. ( p.ll| ) one can get the following solutions 



d 

5ju[s - U{x),x] = 5U^{x)—j{e) + Aj(e). 



(5.^ 



5Ue - U{x),x\ = j^ die5E^{i)—f[e - U{i),i] 



eX[e - urn ^^'^^ 

which show that Sf has the above mentioned property. Here Aj{e) are the fluctuations 
of the current at the left boundary x = 0. The fluctuations of the distribution function 
A/(e) at the right boundary are assumed to be zero. If we assume X{e) to be a constant 
(not depending on the energy), taking into account Eqs. (5^) and ( |5.9|) and the equation 



Sf^{e, 0) = we immediately arrive at the result obtained by Nagaev [|T^ 

AJ = jf dx I deg[e-U{x),x\. (5.10) 



B. Comparison with approach of Ref. [4] 



Now we embark on setting forth the crucial point of the paper. Eq. ( [5.7| ) does not 
coincide with the equation for the field fluctuations obtained in |TT| by direct linearization 
of Eq. (O) for s = 



d 
dx 



5U^{x)^E{x) 
dx 



d 
dx 



U{x)^6EM 
dx 



+e^±E{x)6EUx) 
2 dx 



A-n 



ADnn 



{6J^ - G^). 



(5.11) 



It is necessary to understand the origin of this discrepancy. 

First, we adopt for the time being the scheme of Ref |Tl| and reconsider Eq. ( [4.3|) for 
the current 



J d f°° 

— = —e—- / deu(e)D(e)f(e,x) 

A dx J-u{x) ^ ' V ;j V 7 ; 

+-Doe^E{x) / deiy{e)fie,x). 
2 J-u{x) 

9 



(5.12) 



For the total current (the d.c. current plus fluctuations) the equation reads 



J + 5J -G 

A 

3 



d f°° 
-e— / de p{e)D{e) [f\e, x) + 5f{e, x)] 

dx J-U{x)-5U{x) 



■-! PC 

+^D,e'[Eix) + 5Eix)] 
z J —I 



deu{e)[f{e,x) + 5f{e,x)] 

-U{x)-5U(x) 

Taking into account Eq. ( p.l2| ) we get the following linearized equation 



(5.13) 



3 { r^ 

+ -Doe'E{x) y_ 



SJ — G d r°^ 

-. — = -e-r- / deu{e)D{e)6f{e,x) 

A dx J-u(x) 

deu{e)6f{e, x) + 6U——— / deu{e)f{e, x) \ 
-u(x) dU[x) J-u{x) J 

+-Doe^6E{x)J_^^ deu{e)f{e,x)-e—SUj^^j-^ / _ deiy{e)D{e)f{e,x) 



(5.14) 



-U{x) • . - • . i^j, 6U{x) J^U(x) 

If one linearized the Poisson equation in the spirit of Ref. |jll[ one would see that the term 
in the curly brackets in Eq. ( p.l4| ) would coincide with {n/A'Ke){d5E/dx), so that 



K d6E 



30 ^ 

deu{e)6f{e, x) + SU——- 
u{x) dU[x) 



dev{e)f{e,x). 



(5.15) 



Ane dx J-u(x) • - • • ou(x] J-u{x) 

Simplifying the first, third and fourth terms on the right-hand side of Eq. ( |5.14D with the 
help of Eq. (|5.6| ) and inserting instead of the term in the curly brackets {n / ATie){d5E / dx) 
we arrive at 



SJ-G 

A 



d I j~, rr/ X '^ d6E\ 3^ K d ^_^ ^ 



dx 



An dx 



An dx 



hi d S 


\ dE] 


SU 


U 


Ivr dx SU 


dx 



(5.16) 



One can see that the last term on the right-hand side of this equation coincides with the 
first term on the left-hand side of Eq. (|5.11|) . To avoid confusion note that we believe 
Eq. (5.15) to be also wrong. We have written it here only for the sake of detailed com- 



parison with the approach of Ref. ITT]. We believe that the correct Poisson equation for 
the fluctuation field is 

K dSE 



dei'{e)6f{e,x) 



(5.17) 



Ane dx J-u(x) 

Now we will add in Eq. ( ^.3| ) for the d.c. current the terms that actually vanish as 
they are proportional to the integrals of the distribution function over e with the upper 
limit —U{x) whereas the distribution function f{e,x) = for e < —U{x). The point is 
that when we calculate the fluctuations by replacement U{x) -^ U{x) + 6U{x) they will 
give a nonvanishing result. We have 

J d 

A 



-U{x) d rco 

—e—- I d6i'{e)D{e)f{e,x) — €—- de i'{e)D{e)f{e,x) 

dx Jo dx J-u{x) 



+-Doe'E{x) 



'U{x) 3 roo 

deu(e)f(e,x) + -Doe^E(x) deu(e)f(e,x). (5.18) 

2 J-u{x) 
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Rewriting this equation for the total current we get 
J + 5J-G d r-u{x)-su 



A 



dx Jo 
d /■°° 



dEiyiE)D{e)[fie,x) + 6fie,x)] 
deu{e)D{e)[f{e,x) + 6f{e,x)] 



dx J-U{x)-5U(x) 
Q l—U{x)-5U{x) 

+-D,e'[E{x) + 6E{x)] / deu{e)[f{e,x) + 5f{e,x)] 



< re 

+ -D,e'[E{x) + 5E{x)] 
z J —I 



'U(x)~5U{x) 

Linearizing this equation and using relations like 



deu{e)[f{e,x) + 6f{e,x)]. 



(5.19) 



SU{x) 



6 



de iy{e)D{e)f{e,x) 



-U{x)-5U{x) 



deu{e)D{e)f{e,x) (5.20) 



6U{x) J-U{x) ' ' ■ ' - ■ ' J-U{x) 

we arrive at Eq. ( |5.1| ) that has been derived above. We see that the cancellation of the 
linear in 5U contributions in Eq. (|5.19|) is due to the terms that vanish in the equation for 
the d.c. current but should be taken into account when one considers fluctuations. This 
is why the linearization of Eq. ( [1.5D leads to Eq. ( |5.11| ) that we believe to be wrong as it 
does not take into account all the sources of fluctuation or, in other words, all the terms 
in Eq. ( |5.18| ) containing U{x). 

The solution of Eq. ( [5.7] ) with boundary conditions 



E{x)5E^{x) |^_^o 
U{x)^5E^{x) 



0, 



(5.21) 



dx 



IS 



ADqk 

At: 



5E^{x) = V^'^x) 



x ^0 



C' + ilTJT^Ty/o*''*^" -«("'' 



where C is an integration constant. Requiring a nonfluctuating applied voltage 

/ dx 6E^ = 
Jo 

we get from Eq.( |5.22D that the constant is 



(5.22) 



C 



(ti 



\iiL) 



^)i7=75^r^«'^""-^"«»' 



where 



ip^x) 



diU'l\i). 



(5.23) 



(5.24) 



Now we require at the right boundary 



11 



;^*^"W 



X = L 



(5.25) 



and get 



where 



1 fL 
SJ^ = — dxIl{x)G^{x) 
Z Jo 



dU'(L)U'^/^(L) fL , xtpix) 
Z = L -\ _ , , ,, / ax ■ 



2ip{L) Jo f/'^/2+i(x)' 



(5.26) 



(5.27) 



n(.)^l + ^^\^y:^(^V^r ^(0 



2^{L) Jx '^^U^/^+^iO' 
Then the noise power P is 

P=— [ dxU^ix) < G^{x) >^ . 
Z^ Jo 

Since, according to Eq. ( |5.5| ) 

< O'ix) >^= 2e^A / deu(6)D(e)f(6,x) = 2e^ADoU(x)- -^ 

JO 47re dx'^ 

Finally we arrive at 



P^^j^f'i.u^.m./" 



(5.28) 



(5.29) 



(5.30) 



(5.31) 



AnZ'^ Jo dx'^ 

The potential distribution can be found following method of Ref. [|ll|, i.e. solving Eq.( [4.5|) 
with boundary condition Eq.(|3) at a; = L. Using Eqs. (|]2D, (jO^ ), (|]28D and ( ^311) 
we calculate the suppression factor P/Ppoisson- For physically relevant different values of 
the dimensionality d we get 



-' / -* Poisson * 



0.3188 for d = 3, 
0.4512 for d = 2, 
0.682 for rf = 1. 



(5.32) 



Thus, in this particular case our results differ from those calculated in Refs. [|r^ both 
analytically (that, in our opinion, is of principal importance) and numerically (although in 
this particular case the difference is not great). Naturally, there is essentially no difference 
with the results calculated within an ensemble Monte Carlo scheme in Ref. |1^ . 
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C. Energy dependent scattering time 



Here we calculate the noise power for s = ±1/2 and d = 3. The equation for the 
fluctuations is 



^(5J.-G.)-- 



'uy 



dx 



e^l±^(-Uy±(E5E^). (5.33) 



Introducing the dimensionless potential x given by Eq. ([4.7|) and the fluctuation of the 
field AE 






l/(s+2) 



AB(| 



(5.34) 



one can rewrite Eq. 



as 



^E"Ul-^\^^E--is + ^V-AE- ' («-") 



2; X 

We assume s = —1/2, d = 3 and get 



2/ X 



X 



i+i 



1^1 



(5.35) 



AE" - ^-^AE' - ^AE = x-'"^^^- 
2X X \J\ 



(5.36) 



This equation differs from that derived in Ref. ||TT| while equation for the potential x 
coincides with 



2;^l/2 



x'x" - x'"^" = 1. 



(5.37) 



To calculate the Green's function of Eq. ( |5.36| ) we need the function iIji{x) obeying the 
homogeneous equation 



2 X X 



(5.3^ 



and satisfying the boundary condition '?/'^|x=o = 0. The second function ^2 obeying the 
boundary condition 4'2\x=l = can be expressed through the functions x and ipi 

x"\i) 



^2(2:) = -^1 



'dr^"'«> 



The solution of Eq. ( |5.36| ) can be written using the Green's function 

G{x, x') = ^.^ \Q{x - x')V'i(a;')^2(x) + 6'(2;' - x)'^x{x)'^i{^)\ 
X \^ ) 

as 



(5.39) 



(5.40) 
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AE= f'dx'G{x,x') ^^\^'] J^\ (5.41) 

JO X' [^)\J\ 



Requiring a nonfluctuating applied voltage we get 



SJ = ^ [^ dx^^U(x) (5.42) 

Z Jo x{x) 



where 



U{x)=Mx) f'd^MO+Mx) rd^MO, (5.43) 

Jx Jo 

, dx-^. (5.44) 

x{x) 

Expressing the correlation function < G'^(x) > through x we get for the shot noise power 
reduction factor 



n-oisson Z^ Jo X^' [X) 

We determine potential x following Ref. ||ll| and find numerically ipi from Eq. ( ^.381) . 



Now the functions ip2, H and the constant Z can be found from Eqs. (|5.39|) , ( ^.43| ), and 
(|5.44|) . The reduction factor can be evaluated as 



P/Ppoisson = 0.4257 (5.46) 

which is about 10% larger than the result obtained in Ref. 0]. As indicated by Gonzalez, 
Gons,lez, Mateos, Pardo, Reggiani, Bulashenko, and Rubz flOl, they obtained for s = 
—1/2 as a result of numerical simulation 

P/Ppoisso„ = 0.42 - 0.44. (5.47) 

One can see that this interval is noticeably nearer to the value given by Eq. (|5.46|) than 
the result of Ref. g^. 

In the case s = 1/2 the reduction factor can be evaluated as 

P/Ppoisson = 0.1974. (5.48) 



that is slightly smaller than the result of Ref. |10 
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VI. CONCLUSION 

In summary, we have developed an analytical theory of shot noise in a diffusive con- 
ductor under the space charge limited regime. We find difference between the present 
theory and the approach developed earlier and indicate a possible origin of the difference. 

Several conclusive remarks. The calculated nonequilibrium shot noise power in a non- 
degenerate diffusive semiconductor for two types of physically relevant elastic scattering 
mechanisms turned out to be very close to the ones obtained in numerical simulations 
by the authors of Ref. |TU|. The computed noise suppression factor P/Ppoisson for the 
case of an energy-independent scattering time is also rather close to the analytical results 



obtained earlier by Schomerus et al. [11 1. However, for an energy- dependent scattering 



the numerical difference between our results and the ones of Ref. 11 is considerable. 



Let us clarify once more the point as to why the authors of Ref. |11| arrived at the 



equations that differ from ours. As an example we take the Poisson equation. According 



to Ref. Ill one could write 



/oo 
deiy(e)f(e,x), (6.1) 

-U{x) 

where n, U are the exact tota/ concentration and potential energy, / is the toto/ distribution 
function(the mean value plus the fluctuating part). Linearization of this equation leads to 
equations of Ref. []TT| . The authors of Ref. |TI[] could have argued that since the voltages 



in the reservoirs do not fluctuate and U is set to zero at the left boundary and since the 
total energy S = e + U is and remains positive therefore the total distribution function is 
zero for e < —U. 

Our point is that one cannot justify Eq. (|6.1|) for the total values of these variables 



including the stationary and fluctuating parts. This is readily seen from the fact that the 
fluctuating part of the distribution function itself depends implicitly on the mean value of 
the distribution function through the correlation function. One should bear in mind that 
an equation involving both the mean and the fluctuating quantities must be regarded 
symbolically. Indeed, such an equation is in fact equivalent to two equations: one for 
the mean values and the other for the fluctuating part. Regarded literally it can lead to 
confusion. For instance, analyzing the equation 



n + 6n= _ deiy{e)(f + 5f) 

J-U-5U 



-U-5U 

one could have come to a wrong conclusion that the mean value n depends on such an 



average as 6U6f. 

A few words about the boundary conditions for the potential. The boundary conditions 



used are not applicable within the length Ry = jK,V/ATTen(0) near the electrodes. As the 
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nonequilibrium noise power is a bulk property [mark, for instance, integration over the 
coordinate in Eq. (|5.45| )] this approximation is justified since we assume that the sample's 



length L is much bigger than Ry- 

Being interested in analysis of the fluctuation phenomena in the simplest situation 
of the space-charge limited diffusive conduction regime, we have not taken into account 
the electron-electron collisions. Meanwhile, such collisions can bring about additional 
electron-electron correlation which one should consider treating a more general case. 
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